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Abstract: In this study, the use of taylor series method in the calculation of some means with single auxiliary variable 
developed in a simple random sampling and mean square error unit ratio estimators having certain properties was 
investigated. An application was performed in respect thereof. The study population (mass) included 111 secondary schools 
from 18 districts of Trabzon province. Auxiliary variable (x) was taken as the number of students whereas the main variable 
(y) was taken as the number of teachers. Sample size was calculated as 45 for unit ratio that has certain features. Afterwards, 
theoretically proposed mean and units ratio estimators having certain properties were compared numerically. Random 
sampling was performed using the SPSS 20 program thus giving an equal chance to the units sampled and variability in the 


population was protected. 
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1. Introduction 


Proportional estimator method is not a sampling method 
but a method for finding estimator. In the case where the 
correlation between two variables is positive, information 
from auxilary variables is used. Newly proposed ratio 
estimators are frequently encountered in the literature. Data 
related to auxiliary variable were examined in detail by 
Cochran (1977). Singh-Tailor (2003) estimation is based on 
mean ratio estimation proposed in the case where mass 
correlation coefficient (p) of main and auxiliary variables 
is known. Kadilar and Cingi (2004) proposed new ratio 
estimators by implementing linear regression estimator 
rather than simple estimates and using coefficient of 
variation information (C,,) from auxilary variable. Kadilar 
and Cingi (2006) developed a new mean ratio estimator by 
typing regression estimation itself in Singh and Taylor 
estimator and using coefficient of mass variable from 
coefficient of mass correlation and auxiliary variable rather 
than simple estimator. Naik-Gupta estimator (1996) is the 
unit ratio estimator having a certain property. Singh et al. 
(2008) proposed an estimation method based on the 


implementation of linear regression estimation having 
certain features rather than simple conventional estimation 
in Naik-Gupta estimation and the use of data such as 
variation coefficient and correlation coefficient. 

In our application, the estimators proposed in the 
literature were compared numerically. 


2. Taylor Series Approach 


Observation of non-linear estimates is frequently 
encountered in sampling. It may not be possible to calculate 
the variance of these non-linear estimators exactly. 
Estimations such as ratio estimate or estimation of 
regression coefficient may be used. To calculate the 
variance of a non-linear estimator, observations of non- 
linear estimator are approximated to a linear function and 
an appropriate variance equation is written. This 
linearization is usually performed by Taylor Series Method. 
In other words, Taylor series method can be used in the 
calculation of variance of an estimator. 
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In general, it can be given for the variable k as follows: 
A(%, Yor Pe) = AM Yo, 


where A(%, Yoo, Yx) is a non-linear function and h(%,, Y,," 
--V#is the estimation of this function. 


Pee Oh(@4,42,"",@K) 
J —_ da; (2) 


Yoni Pah yp) rk 
where 0 represents the notation of partial derivative and 


1 07A(%4,Yo,-~ 


Yi rg oS 
Ric(Vir @) = Tihs Dla 5 ee (G — YAK) + Oe 


(3) 


Here, O, represents the terms after second order in Taylor 
series method. 

Here, the term R;,(Y¥,, a) does not change the value of 
mean square error dramatically; but R,,(¥,, a) is important 
to calculate the bias. In simple random sampling, second 
order terms, in other words, the term R,(Y¥,,a@) is 
frequently neglected when the sample size is large enough. 
Therefore, the mean square error can be obtained by Taylor 
Series Method as follows; 

A(Y, Yor Me) & AM Yor Me) + UL a (% -Y) A 

By using Equation (4), mean square error is obtained as 
follows; 


E[h(Py, Py, %) — AM Yo,.Y]) = E[zvi1d(% - ¥)| (5) 


MSE[h(h, Yo FO] = ELad%-Yyl = © 


3. Mean Square Error of the Proposed 
Estimators 


3.1. Simple Conventional Estimation and Mean Square 
Error 


When n-width sample is taken from n-width mass via 
simple random sampling method, variance of ¥y can be 
obtained as follows; Variance of mean estimator, 


2 
(Won) Sy? 
N n 


VQ) =EG-Y) = (7) 


Where f = ve A= is taken, Equation (8) becomes 


_, _ (N-n)Sy? _ 1-f _ 2 
Vy) Se ASy (8) 
22 Diedink? 
a= yee 
3.2. Conventional Propotional Estimation and Mean 
Square Error 


Conventional ratio estimation for (Y) mean mass of 
variable of interest y is 


Ye) + Dh di (% —¥)+ReMe@) (1) 
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Yo = 3X = RX (9) 


Here it is equal to R = and y is the sample mean of the 


variable of interest and xX is the sample mean of the 
auxiliary variable. We assume that mass mean of auxiliary 
variable x is known, 


x= Det 


where n is the number of sample units. 
Mean square error of this proposed estimator was 
obtained as 


xX, VEY = ya1 i (10) 


MSE (J) © A(Sy? — 2RpyxSySq + R2S,7) = A(S,? + R?S,2(1 — 20)) (11) 


Where = pz, C, and C, are the variable of interest and 


coefficient of variation of the auxiliary variable respectively. 


A= a is the sampling fraction f and it is obtained by 


division of sample width to the mass width. In other words, 


n 3 - ¥ 
os the number of N mass units. R = z is the mass 


ratio; Ss and S,’are the mass variance of the variable of 
interest and auxiliary variable respectively [2]. 


3.3. Investigation of Singh and Tailor Ratio Estimation 
and Mean Square Error 


If the mass correlation coefficient of variables of interest 
and auxiliary variable (p) is known, ratio estimate 
proposed by Singh-Tailor can be defined as follows; 


7. = 2_(X 
Ver =o TP) (12) 
Mean square error of this estimator is expressed by 
MSE (Fer) = A (sy + R25,2w(w — 26)) (13) 


xX pc. 
where w = =— ve 0 = —[3]. 
X+p Cx 


3.4. Some Estimators Proposed by Kadilar-Cingi and 
Mean Square Error 


Kadilar and Cingi (2004) proposed the following ratio 
estimators. 


(14) 


Y+b(X-X) 
Tae, = SX 


= ¥+b(X—x) 
Fos =, K+Gl (15) 


Riel = Toe x) 
variance of acting variable and sy, is the sample 
covariance between variables of interest y and auxiliary 
variable x. General form of this expression was proposed 
by Ray and Singh (1981) [4]. Proportional estimator 
proposed by Ray and Singh is as follows; 


Ss 
and b= =o Sy” represents the sample 
x 


isa xe 


Yrs = (16) 
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In this estimator, if b=0 and §=1, it becomes 
conventional proportional estimation. In addition, in the 
case of a = 1 and § = 1, proportional estimator proposed 
by Ray and Singh corresponds to that proposed by Kadilar 
and Cing1. Mean square error of Kadilar and Cuingi 
estimator is expressed as follows; 


MSE(xci) = AL(B + Rxci) Sx? + Sy? — 2(B + Rxci)Sxy)] (1 


7) 
Here, if B = aye and is squared 
Sx 

It is expressed as 

MSE(Fxci) = AL(Reci? Sx? + Sy? — p?)] i = 1,3 (18) 
Y 

where Ryq = R = rae Rx = me, 1 

Inference: In Equation (17), if B = —Rxcj, mean square 


error Of Yx¢j estimator is equal to the y variance of simple 
estimator. 
If B= —Rxci, 


MSE(Fxci) = AS,” = V(¥) 


Kadilar and Cingi (2006) developed a new ratio 
estimator by utilizing regression estimation rather than 
simple estimation in Singh and Taylor Estimator and using 
coefficient of mass correlation and variation coefficient of 
auxiliary variable. This estimator is as follows; 


= V+b(X-*X) -s 
Fors = +p) (19) 
a V+b(X-X) -s 
Vora = a AC, +p) (20) 
= V+b(X-X) zs 
Vora = Fong, Ap + Cx) (21) 


Mean square error of this estimator can be expressed as 
follows 


MSE (Ypri) = A[Rpri’ Sx? + Sy?(1 — p?)] i= 1,2,3 (2 


[6]. 


— __¥ _ _¥p 
where: Hag X+p Rpra = XCytp’ ~ PT3 Xp 4+cy 
3.5. Some Ratio Estimators that have Certain Feature 

which were Developed in Simple Random Sampling 
Using Auxiliary Information and Mean Square Error 


Here, R in ratio estimation is replaced by B in regression 
estimation. The use of auxiliary variables can improve the 
efficiency of estimation in the case of where auxiliary 
variables and the variable of interest are highly correlated. 
Suppose that n-width sample is taken from n-width mass 
(population) via simple random sampling without 
replacement. y; and d; represent the observation values in 
the variables for i = 1,2, ..., N respectively. Suppose that 


1,if the desire feauture is achieved 
oa {rt , (23) 


0, other cases 


: , ; A 
By using this Equation, P = rc and p =< mass and 
sample ratio estimations and ratio equation of the desired 
events can be obtained. Here, = Y%,; and a= Yo} 
show the total number of mass and sampling units 
respectively [1]. 


3.5.1. Investigation of Naik and Gupta ratio Estimators 
and Mean Square Error 


Tne =>P (24) 


where y is the sample mean of the variable of interest. p 
is the unit ratio having certain features in the sampling. 
Mean square error of this estimation is 


MSE (Ye) = A(Rne’ So” — 2RyoSyg + Sy") 


where Ss = rep — P)* and 
1 —_ 
Syg = = dizi — P)G% — Y) (71. 


3.5.2. The Proposed Singh, Chaucer, Sawan and 
Smarandache Estimates and Mean Square Error 
Error 

Estimations following Ray and Singh Ratio Estimations 
are 


_ ¥+by(P-P) 


aloud (25) 
_ ytbg(P-p) 

ae fem [P+C,| (26) 

ty = OE TP + pol (27) 


PtPo 


In Equation (25), if by =0, the proposed estimator 
becomes Naik and Gupta estimation in Equation (24). 
Mean square error of this estimator is 


MSE(t;) = A[Ri?Sg7 + Sy7(1 — p07)] 5 i = 1,2,3 


4 4 
R3 


Y 
Ry =e R2 = > = 
P P+Cp P+po0 


where and 


So? = Tb: — PY 
Syp = —— Dh: — P)(y; — Y) [8]. 


4. Theoretical Comparisons 


4.1. Comparison of Singh and Tailor Estimator and 
Conventional Proportional Estimation 


Conventional proportional estimator; 


MSE (Jo) = A(R?S,? — 2RpyxSySq + Sy") = A(S,? + R?S,2(1—28)) (28) 


We have seen conventional proportional estimator in 
Equation (28). Conventional proportional estimation is 
obtained by inequality of mean square error of Singh and 
Tailor Estimator 
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A(Sy? + R?S,’w(w — 20)) < A(S,? + R?S,7(1 — 26)) 
w(w — 20) < 1-20 
w?-—1-2w0+20<0 
(w —1)(w +1) —-20(w-1) <0 
(w-—1)(w+1-26@) <0 

This inequality is satisfied when 
w-1<0 vew+1-26>0 (29) 
or 
w-1>0vew+1-20<0 (30) 


Under these conditions, mean square error of Singh and 
Taylor Estimator is less than that of conventional 
proportional estimator. Therefore, it is more effective. 


4.2. Comparison of Singh and Tailor Estimator and 
Simple Conventional Estimation 


a (ee + R2S,2w(w — 26)) ZAS,° 
R?S,.?w? < R?S,?w20 
Ww 
awe (31) 


Under condition (31) Singh-Tailor estimation is more 
effective than simple conventional estimation. 


4.3. Comparison of Kadilar-Cingi (2004) Estimation and 
Conventional Proportional Estimation 


From Equation (11) and (18); Under the condition 
A[(Reci?Sx” + Sy?(1 — p?)] < A(Sy? — 2RpSyS, + R7S,”) 
(Rui? — R?)5,7 < p?Sy? — 2RpSySz) 
(Reci’ — R?)Sx" < pSy(pSy — 2RS,) (32) 
Kadilar-Cing1 estimation is more effective. 


For example, consider the case i = 1 


Y 
If Rygi = R= 5 


Under the condition 


0 < pS, (pS, — 2RS,) 2RSy < pSy 


peo (33) 


Cy 
Kadilar-Cing1 estimation is more effective. 


4.4. Comparison of the Proposed Singh, Chaucer, Sawan, 
Smarandache Estimations and Simple Estimate 


Comparison of the proposed estimators and simple 
estimates and variance of the sample mean was examined 
in simple random sampling as follows; 


V(¥) =AS,? (34) 
V(y)-—MSE(t,))>0 ;i=1,2,3 (35) 
ASy? > A[R?S47 + Sy?(1 — p07)] (36) 
2. 56° p2 
Po > 5 2Ri (37) 
y 


condition is obtained. Under this condition, the proposed 
estimator is more efficient than simple estimates [8]. 


4.5. Comparison of the Proposed Estimators with Each 
other 


In the case of 
A[Ry? Sy? + Sy?(1 = po?) < A[R;7*Sg" + Sy? = po), fF =2,3 (38) 
RES Ry oa 
R,?S¢" — Ry’ Sg" <0 
So’ (Ri? — RB?) <0 
IRs| < |R,| (39) 
t, estimate is more effective than t, and ts. 


4.6. Comparison of the Proposed Singh, Chaucer, Sawan, 
Smarandache Estimates and Naik and Gupta 
Estimations 


With the proposed estimation, comparison of Naik and 
Gupta estimate is achieved under the condition of 


MSE(jyc) — MSE(t) >0 3i=1,2,3 (40) 


Sq? Cc 
p> sF Ri? — Rg’ + 2RgPo = (41) 


Under this condition, mean square error of the proposed 
estimator is less and therefore more effective [8]. 


4.7. Comparison of Naik and Gupta Estimate and Simple 


Estimation 
ASy* > A(Rug? Sg” — 2RucSyo + Sy”) (42) 
Dp SySe°>' Rueda (43) 
RnGS@ 
>t 44 
Po > a5, (44) 


condition is obtained. 

Under this condition, we can conclude that mean square 
error of Naik and Gupta estimation is smaller and therefore 
more effective. 


5. Application 


School ranking of 18 districts in Trabzon province in the 
2011-2012 academic years was used (Trabzon MEM, 2011- 
2012) [9]. The data from 18 districts of Trabzon province 
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included the number of students and teachers from The relationship between the number of teachers and the 
secondary schools in the 2011-2012 academic years. number of students can be demonstrated by a linear 
Information regarding the number of teachers and students equation (y = bx) passing through the point of origin and 
was obtained from Trabzon provincial Directorate of the related variables can be estimated using auxiliary 
National Education. variables. Mass unit is the 111 secondary schools in 

Main variable (y) is taken as the number of teachers and Trabzon province. Data regarding the mass is given in 
auxiliary variable (x) is taken as the number of students. Table 5.1. 


Table 5.1. Mass units regarding the number of teachers (vy) and the number of students (x) variables. 


N=111 Y = 29.28 X = 395.06 p = 0.939 
pene S,? = 651.312 Sy = 25.521 C Dae 
os i Yy y y Y 5 

2 fe Sx Sz 
S,2 = 159948.937 S, = 399.936 Cy = == 1.012 Cyn = Fp = 0.829 


Samples are obtained from a mass via simple random of the schools in the city. where d is the amount of 
sampling without replacement using SPSS-20 program. tolerance. Considering the above equation (45), sample size 


First, sample size is calculated. reduces when the the amount of tolerance increases. The 
’ amount of tolerance (d) depends on the researchers. In the 
ie (45) present application, if d = 0.09 sample size is obtained 


~ d2(N-1)+t2PQ - 


can be used in the estimation of sample size [10]. 

In this equation, P ratio of units having the same features 
is (Q = 1-—P). If we do not have information about units 
ratio which has the same feature, the largest sample that can The statistics and the data regarding the schools location 
be selected for P= Q=0.5 could be reached. However, _ in the city are given in Table 5.2 and Table 5.3 respectively. 
unit ratio having specific features is calculated as the ratio 


i 1,967*111+0.81*0.19 
~ 0,092*110+41.962+0.81+0.19 


45 (46) 


Table 5.2. Sample units of the number of teachers (y) and the number of students (x) variable. 


n=45 y = 32.02 x = 413.82 r=0.955 
Ry n 
k= La = 0077 See = 875.522 Sy = 29.589 f= N = 0.405 
DG 
2 eth Se _. 
Sy, = 183282.5586 S, = 428.115 A=— = 0.013 b= eps 0.066 
n x 


Table 5.3. Mass and sample data required for unit ratio having certain features. 


5 
A=90 a=39 P =0.810 Po = =e = 0.216 
yr 
by = 22% = 14.006 S,? = 651.312 S, = 25.521 Sy 
pr gen is 2 = 651. eee: Cy = 7 = 0.872 
Ss S 
Sp? = 0.155 Sy = 0.394 Cy = a = 0.486 Cn = sd = 0.092 


Table 5.4. Mean Ratio Estimators and Mean Square Error 


Mean Estimators Mean Values Mean Square Error 
Simple Estimation (7) 32.02 8.467 

Conventional Proposional Estimation (y,) 30.568 2.168 

Singh-Tailor Estimation (57) 30.572 2.163 
Kadilar-Cingi-1 Estimation (¥xc1) 29.386 11.190 
Kadilar-Cingi-3 Estimation (Vxc3) 29.389 12.388 
Kadilar-Cingi Estimation (Vp,1) 29.389 12.357 
Kadilar-Cing Estimation (Vp,2) 29.389 12.082 
Kadilar-Cing Estimation (Vp,3) 29.390 13.880 


According to the mean ratio estimators, the average province was 30. In addition, Singh- Tailor estimation 
number of teachers per secondary schools in Trabzon (Ys) showed the highest sensitivity. 
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Table 5.5. Ratio Estimators having certain features and Mean Square Error. 


Estimators of units ratio having certain features Estimates Mean square error 
Naik-Gupta Estimate (Vy) 29.949 9.058 
Singh-Chaucan-Sawan-Smarandache 1| Estimate(t, ) 29.216 10.705 
Singh-Chaucan-Sawan-Smarandache 2 Estimate(t,) 29.940 9.10 
Singh-Chaucan-Sawan-Smarandache 3 Estimate(t;) 29.619 9.713 
According to ratio estimators having certain units feature, For this condition i = 1; 
the average number of teacher per secondary schools in 
Trabzon was 29. Naik-Gupta Estimate (jy,) estimator was 2 2 Sp”, 2 _ 0-155 * (36.148)? 
aan © One Po? = 0.216? = 0.046 <  R,? = ——_—_——_ = 0.311 
the most efficient estimator with the minimum mean square Sy 651.312 


error. 
Comparison of Singh and Tailor Estimataion with 
Conventional Proportional Estimation 
MSE (Ys) = A(Sy” + R?S,’w(w — 28)) = 2.163 
where w = = 0,998 and @ =" = 0.809 
X+p C 


To ensure inequality MSE(¥s7) < MSE(¥o) w—1 <0 
and w+1-20>0 orw-—1>0andw+1-20<0 
condition should be fulfilled. Thus, since w — 1 = 0.998 — 
1<0O and wt+1-20=0.998+ 1-2 * 0.809 = 
1.967 > 0 condition is satisfied, 


MSE(¥sq) = 2.163 < MSE(¥o) = 2.168 


IR 


mean square error of Singh and Tailor Estimation is smaller 
and therefore more effective. 

Comparison of Singh and Tailor Estimation and Simple 
Estimation 

The variance of simple estimation is V(y) = AS? = 
8.467 

MSE(¥sr) < MSE(¥) inequality is valid when the 


condition i —— = 1.234 < 2 is satisfied. Since this 
condition is fulfilled, the Tailor Singh estimation is more 


effective than simple classical estimation. 


Comparison of Kadilar-Cingi-1_ Estimation and 
Proportional Estimation 
Under the condition = MSE(¥c1) = A[(R?S,7 + 


Sy21-p2=11.190 p>2CxCy, Kadilar 


Cingil estimation is more effective than conventional 
proportional estimate. However, condition is not satisfied, 
so conventional proportional estimate is more efficient. 

Comparison of Naik-Gupta Estimation and Simple 
Estimation 


Mean square error is found as MSE(jyg) = 
A(Rya? So” — 2RvcSyg + Sy?) = 9.058. When compared 
with simple estimates, if the condition p, > “Nee is 

y 


satisfied, Naik and Gupta Estimation is smaller and 


therefore more effective than simple estimation. However, 
; RncS are : 
since Py = 0.216 < — condition is not satisfied, error 
y 


of simple estimate is smaller and therefore more effective. 
Comparison of the Proposed t, Estimation and Simple 
Estimation 


To ensure V(y) — MSE(t;) = 0 
2 


“6; Ri? condition must be fulfilled. 
'y 


inequality, p,* > 


; Sy? Se : 
Since p,* > Ry? ondition is not satified, error of 
y 


simple estimation is smaller and therefore more effective. 
For i = 2; 


2 
po? = 0.2162 = 0.046 < =f Ro? 


0.155+*(22.593)? 
gia, ae 
0121’dir. Therefore, simple estimation is more effective. 
In the same way, simple estimate is more effective for 
i = 3 case. 

Comparison of the proposed t, estimation and ty 
estimation 


Since Ry =~ = 36.15 > R, == = 22.59, 
Pp 


MSE(t,) = 9.10 < MSE(t,) =10.71 error of t, 
estimation is smaller and therefore more sensitive. 


6. Conclusion 


When B= —Rxqj in Kadilar-Cing: (2004) and B= 
—Ry,; in Kadilar-Cingi (2006) estimation, mean square 
error of Kadilar-Cingi (2004) and Kadilar-Cingi (2006) 
estimates are equal to that of simple estimation. 

When comparing Singh and Tailor Estimator with the 
conventional proportional estimators, under the condition 
w-1<0 and w+1-28>0 or w-1>0 and 
w +1-— 26 <0 Singh-Tailor estimation is more effective 
than conventional proportional estimation. 

In our numerical example, since 
w—-1=0.998-1<0 and w+1-26=0.998 + 1— 
2 * 0.809 = 1.967 > 0 condition is satisfied, MSE (V5) < 
MSE (Yo) inequality is obtained. 

In the comparison of Singh and Tailor Estimators with 
the simple estimate, aS 2condition is satisfied and under 


this condition, Singh-Tailor Estimator is more effective 


Tae w 0.998 
than simple estimator. In our applications, > = 5 = 


1.234 < 2 condition is fulfilled and Singh-Tailor Estimator 
was found to be more sensitive. 

In the ratio of units having certain features, estimators 
proposed by Singh et al., t,, t2 and t3 estimators were 
compared and R, < R; (j = 2,3) condition was obtained. 
Under this condition, the proposed t, estimation is more 


effective than t, and tz estimations. In our application, 
4 Y F : 3 
R, = a 36.15 > Ry = Pep = 22.59 inequality is 
valid and expressed as MSE(t2) < MSE(t,). 


Considering the numerical examples, according to the 
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estimators proposed in the literature, the average number of 
teachers per secondary schools in Trabzon province in 
2011-2012 academic years ranged from 29 to 32. Singh- 
Tailor estimation (sr), is the mean ratio estimator having 
the smallest root mean square error (2,163) whereas the 
estimator proposed by Kadilar and Cingi (¥p-3) has the 
greatest mean square error (13,880). Among the estimators 
having certain features, Naik-Gupta estimator (9.058) has 
the smallest mean square error. 

The numerical results obtained are valid only for this 
application. 
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